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Abstract. We present here a detailed study of the behaviour of a three dimensional Brownian motor 
based on cold atoms in a double optical lattice [P. Sjolund et a/., Phys. Rev. Lett. 96, 190602 (2006)]. 
This includes both experiments and numerical simulations of a Brownian particle. The potentials used are 
spatially and temporally symmetric, but combined spatiotemporal symmetry is broken by phase shifts and 
asymmetric transfer rates between potentials. The diffusion of atoms in the optical lattices is rectified and 
controlled both in direction and speed along three dimensions. We explore a large range of experimental 
parameters, where irradiances and detunings of the optical lattice lights are varied within the dissipative 
regime. Induced drift velocities in the order of one atomic recoil velocity have been achieved. 

PACS. 32.80.Lg Mechanical effects of light on atoms, molecules, and ions - 05.40.Jc Brownian motion - 
32.80.Pj Optical cooling of atoms; trapping 



1 Introduction 

Brownian motors (BMs) are small scale engines that con- 
vert random fluctuations into deterministic work. In or- 
der to realise such a device, symmetry has somehow to 
be broken (Curie's principle pQ), and the system has to 
be out of thermodynamical equilibrium [2]. Much effort 
has been invested in studying the underlying mechanisms 
of BMs operating with the ratchet effect [2,3,4,5 . The 
ratchet mechanism consists in breaking the spatial and/or 
the temporal inversion symmetry of the system so that 
directed transport emerges, using fluctuations as the rele- 
vant input. The paradigmatic device is Feynman's famous 
ratchet and pawl machine [6], based on an idea of von 
Smoluchowski [7]. 

During the past years there has been an extensive body 
of work on BMs and ratchets, several realised in systems 
of periodic arrays of ultracold atoms, known as optical 
lattices (OLs) [8,9, 10, 11( 112]. Optical lattices are periodic 
arrays of micro-traps created by the interference between 
two or more laser beams p~3j[14]. It has been shown that 
ultracold atoms stored in OLs can be controlled and ma- 
nipulated with a very high degree of precision and flexi- 
bility and they are routinely used for example in studies 
of Bose-Einstein condensation and applicable for quantum 
state manipulation [T5([T6]. 
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In this paper, we present an experimental and numer- 
ical study of fluctuations rectified into directed motion 
of Brownian particles in three dimensional OLs [9 j. This 
works by purely optical fields, where ultracold atoms are 
optically pumped between two state-dependent and spa- 
tially overlapped OLs [T7( [T8], both spatially and tempo- 
rally symmetric. These OLs are coupled by optical pump- 
ing, with strongly asymmetric transfer rates, via the vac- 
uum field reservoir, which, together with a shifted relative 
spatial phase between the OLs, causes atoms to be pro- 
pelled in a controllable direction. 



A general understanding of the dynamics of BMs is of 
fundamental interest. In biology, directed transport and 
molecular motors are both driven by stochastic motion, 
thus providing the energy input for the ability of living 
cells to generate motion and forces, e.g., for mobility, con- 
traction of muscles or material transport, and are in them- 
selves a fundamental mechanism of the origin of living 
cells. For example, biological motor proteins which move 
along linear filaments can be described by stochastic mod- 
els coupled to chemical reactions [19 j. So-called ratchet 
models further explain the generation of directed motion 
on the microscopic level out of Brownian motion. The gen- 
eral concept of our BM, based on the idea of [20], may be 
transferable into fields such as chemistry and biology as 
a tool for studies and control of, e.g., molecular motors, 
investigations of intra-cell motion and possible studies of 
Brownian motion in biological membranes [2Tj . 
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Fig. 1. (Colour online). Rectification mechanism. Atoms move 
in two symmetric potentials Ua(z) and Ub(z) that are cou- 
pled with asymmetric optical pumping rates 7a^b and 7b ^ a 
(7a^b ^> 7b^a)- (a) The potentials are in phase. The trans- 
fer from the long-lived state B, to the transient state A, and 
back, does not lead to biased motion, (b) A relative phase shift 
(p is introduced. Spatial diffusion is strongly facilitated in one 
direction, leading to a drift velocity Vd z • (c) The potentials are 
out of phase, ip z = 7r, and no biased motion is introduced. 

2 Concept 

The interaction between atoms and the interference pat- 
tern from the laser beams creates spatially periodic poten- 
tials (micro traps) seen by the atoms due to a second-order 
interaction between the atomic dipole and the light field 
[T3|[14] . Tuning the light frequency of an optical lattice 
close to an atomic resonance provides a dissipative chan- 
nel resulting from spontaneous emission. The OL is then 
accompanied by an efficient cooling mechanism p~4j[22j[23l 
[24] , providing friction in the system, and the presence of 
dissipation, albeit small, will result in a slow normal dif- 
fusion of the atoms in the lattice [T4[ l25 l l26] . 

The basic idea and underlying principle of our BM is 
quite general [20 . Consider a Brownian particle station- 
ary for a short time in one of two spatially and temporally 
symmetric sinusoidal potentials Ua(z) and Ub(z), see fig- 
ure [I] and having a low kinetic energy compared to the 
modulation depths of the potentials. This particle will un- 
dergo Hamiltonian motion around the potential minimum. 
The motion is interrupted by a dissipative process, result- 
ing in a small overall Brownian motion, where the particle 
can be optically pumped between the two potentials. One 
of the components in the asymmetry that eventually gives 
rise to rectification is caused by unequal transfer rates be- 
tween the potentials. Also, we assume the energy damping 
to be much smaller than the typical oscillation frequency 
in the potential wells, so that the dynamics falls into the 
low-damping regime, since no induced drift can be ob- 
tained for sinusoidal potentials in the overdamped regime 

E3- 



By letting the relative spatial phase between the two 
potentials be slightly shifted, combined with the pronounced 
difference in the transfer rates, the spatiotemporal sym- 
metry of the system is broken and a directed motion is 
obtained. The concept can therefore be adapted to induce 
directed motion in any direction in three dimensions since 
the same physical process will function in any direction in 
a three-dimensional case. Because we use symmetric po- 
tentials, the direction of the induced drift can be reversed 
by a proper choice of the relative spatial phase shift. This 
is to be contrasted with the usual use of a ratchet poten- 
tial (e.g., using a sawtooth potential), where the direction 
is fixed by the shape of the potential. 



3 Experiment 

3.1 Experimental setup 

The BM is realised, using cold caesium atoms in a double 
optical lattice (DOL) [17,18,28 j. A magneto-optical trap 
(MOT) is loaded from a chirped-slowed atomic beam, pro- 
duced in a thermal source which fills the MOT with ap- 
proximately 10 8 atoms at a peak density of 10 11 cm -3 . 
Approximately 80% of the initial number of atoms from 
an optical molasses of a few microkelvin is loaded into the 
DOL with a filling fraction of about 0.05 atom per site. 
After loading, the OL light is left on for a chosen time r. 
To avoid spurious spatial drifts of the atomic cloud caused 
by light pressure, great care has been taken to balance the 
laser beam power [28]. Also, it is crucial that the relative 
spatial phase is kept constant. 

This three-dimensional DOL is composed of two spa- 
tially overlapped OLs with identical topography, but which 
can be controlled individually in terms of both well depths 
and relative spatial phase. Two different beams of frequen- 
cies c^a and ub are split into four branches with equal 
power. This lattice geometry is a generalisation of the ID 
lin±lin configuration to 3D [14] , see figure [2^. We use fre- 
quencies that are near-resonant with the D2 line of cae- 
sium at 852 nm reaching the excited fine structure level 
6p 2 P3/2- Each optical lattice operates from a different 
hyperfine ground state of the 6s 2 Si/2 level. One operates 
within the closed F g = 4 — > F e = 5 transition, trapping 
atoms in the F g = 4 ground state. The second lattice 
traps atoms in the F g = 3 ground state, operating in the 
open F g = 3 — > F e = 4 transition, see figure [2Jd. These 
resonances are much narrower than the difference in the 
laser frequencies, which allows us to address the two lat- 
tices independently. At the same time, the wavelengths 
are close enough, and the trapped atomic cloud is small 
enough (typically 1—2 mm in diameter), to ensure that the 
periodicity is the same within the sample volume (e.g., it 
takes about 3.3 cm for the lattices to phase out by tt in the 
horizontal x and ^/-directions). Moreover, we can individu- 
ally control the optical pumping rates between the lattices 
to some extent [18] by tuning the irradiances Ia-, Ib and 
the detunings Z\a, of the lattice beams, which in turn 
changes the well depths, optical pumping rates, diffusion 
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Fig. 2. a) Double optical lattice beam geometry, realised by 
four branches, where each branch contains two spatially over- 
lapped laser beams A and B. ei to are the polarisation vec- 
tors perpendicular to the plane of propagation, b) Schematic 
energy level diagram: Two ground states F g = 4 and F g = 3 
are connected to the excited F e = 5 and F e = 4 states respec- 
tively by two laser fields B and A detuned by Ab and Aa- 
r/27i — 5.22 MHz is the natural linewidth of the transitions. 



and friction, and consequently, also the behaviour of our 
BM. 

By letting the two lattices be out of phase, e.g., in 
the vertical z-direction, while simultaneously having an 
asymmetric transfer rate between the lattices, the motion 
of atoms is channelled in that particular direction [9] . This 
process can be controlled in all three dimensions, with 
regard to both speed and direction, by a proper choice of 
the relative spatial phases ip x , tp y , (p z , irradiances I a and 
2b and the detunings A a and Ab of the laser light. 

The control of the relative spatial phase shifts between 
the lattices in three-dimensions is achieved by extending 
or shortening the lattice beam branches 1, 2, 3 and 
see figure [2^. A relative phase shift only in the horizontal 
^-direction corresponds to an extension of branch 3 and a 
shortening of branch 4 by the same amount, or vice versa. 
In y, an extension of branch 1 and shortening of branch 2 
by the same amount, or the other way around, will cause 
a shift. For a phase shift only in the vertical z-direction, 
a simultaneous extension or shortening of both branches 
1 and 2 or branches 3 and 4 by the same amount will 
work. Combining these, an arbitrary combination of <p x , 
cp y and (p z can be set. The relative spatial phases between 
the lattices are experimentally determined by the kinetic 
temperature, measured by a ballistic time-of-flight (TOF) 
technique [29 , when changing the relative spatial phase, 
see figure [3] A Gaussian fit to the TOF-signal also provides 
the peak arrival time of the distribution and the number 
of atoms by the area under the fitted curve. 



1 This is done by manually adjusting the distance between 
two facing prisms, mounted on linear translation stages in each 
of the four branches. 
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Fig. 3. Measured temperature, T, dependence as a function of 
the beam path extension 5L, which determines the conversion 
between the relative spatial phase to beam path extension [18] . 

3.2 Calibration of the relative spatial phase by 
temperature measurement 

The temperature dependence is shown in figure [3] for a 
relative spatial phase shift in the vertical z-direction. A 

3.3 cm extension or shortening of the beam branches cor- 
responds to a relative spatial phase of 2tt. This is repeated 
for the other directions in order to find the origin of all 
the relative spatial phases. An simplified explanation to 
the temperature variations [17 is an increase in the num- 
ber of scattered photons when a a + point in lattice A 
overlaps with a a~ point in lattice B. Figure^, shows a 
semi-classical model where a <r + point in lattice A over- 
laps with a <j + point in lattice B. If an atom is optically 
pumped from lattice B to lattice A, it will end up in the 
lowest light shift potential (due to conservation of angular 
momentum), and then return to lattice B. In this process, 
an atom only scatters a few photons. The opposite situa- 
tion, where a <j + point overlaps with a a~ point is shown 
in figure [4}d. Here, an atom in lattice B that is pumped to 
lattice A will end up in the least light shifted anti-trapping 
potential. The atom will slide down the potential, gaining 
kinetic energy and, as it is close to a a~ point it will scatter 
photons, which pumps the atom to the lowest potential. 
It is then pumped back to lattice B in a similar process. 
Thus, in this case the number of scattered photons is con- 
siderable higher, which result in heating [T71H8]. In a more 
realistic picture, there is a number of Zeeman sub-levels, 
all corresponding to different light shift potentials. For 
a a + ja~ site overlap, photon scattering increases since 
atoms are optically pumped between a range of different 
M-states. 

3.3 Measurement of drift velocity by TOF 

In figure [5j three TOF signals are shown, all for r = 350 
ms. The middle one, indicated by (0), is when the two OLs 
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Fig. 4. a) If two cr + points in lattice A (ground state |gA)) and 
lattice B (ground state |gB)) overlap, atoms can be transferred 
between the lattices with minimised light scattering, b) In the 
opposite case, when cr + overlaps with a - , heating effects are 
enhanced. 




Fig. 5. (Colour online) Normalised TOF-signals with arrival 
time t as the horizontal axis, for relative spatial phases that 
generates maximum induced drifts in the upward (+z) and 
downward (— z)-direction. The middle one (0) shows a TOF 
signal for no induced drifts when the optical lattices are in 
phase. All three TOF signals are normalised to make compar- 
ison easier and the signals are averages of four drops. 



are in phase and no induced drift is present. Those indi- 
cated by and (—z) are when ip z is set to generate a 
maximum drift, either upwards or downwards. A more di- 
rect technique, and a complement to the TOF technique, is 
imaging the shadow of the atoms, transiently illuminated 
by a weak resonant probe beam, on a CCD detector, see 
section 



3.4 The images provides the centre-of-mass posi- 



tion of the atomic cloud in the xz plane, determined by 
Gaussian fits to the images after an arbitrary lattice time 

T. 

In order to experimentally investigate the behaviour 
of the BM, we change the parameters A a and Ab-, cover- 
ing the frequency span between the excited F e — 3 to the 
F e = 4 states, and the F e = 4 and the F e = 5 states re- 
spectively. Ia and Ib are set independently up to about 20 
mW/cm 2 . To characterise the drift velocity in the verti- 
cal z-direction (i>d z ), we tune the parameters (A a, ^b, ^a? 
Ib) in a number of steps, all for a fixed lattice time r of 350 
ms. For each set of parameters, the relative spatial phase 
in the vertical z-direction (<p z ) was incrementally changed 




Fig. 6. Induced drift velocities Vd z as a function of relative 
spatial phase ip z for three different sets of parameters. Filled 
circles: Ab = 40i^, filled squares: Ab = 25.T and circles: Ab = 
10r, for a lattice time r of 350 ms, A B = 14T, I A = 0.8 
mW/cm 2 and Ib — 1.2 mW/cm 2 for all cases. 



in about 50 steps, covering slightly more than 27r, while 
(f x and cp y were kept at zero, meaning that only vertically 
induced drift was present, either upwards or downwards. 
In total, we measured about 14000 velocity distributions, 
due to the large parameter space. For each change of tp z , 
five TOF-signals were averaged to increase the signal-to- 
noise ratio and Gaussian fits were performed to read out 
the peak arrival time of the TOF signals. 

With a constant drift velocity [9 , the vertical drift 
velocity v^ z is determined from the peak arrival time t of 
the atoms at the TOF probe by 



gt 2 - 21 
2(* + r)' 



(1) 



where I is the vertical distance that the atoms fall down 
to the probe and g is the gravitational acceleration. In fig- 
ure [6| three different curves are plotted, showing v& z as a 
function of tp z for three different values of Z\b, while A a, 
I a and Ib are kept fixed. Here, a clear variation of v^ z 
is evidenced for different values of Ab . In figure [7| Ab is 
plotted as a function of v& z for a relative spatial phase that 
generates a maximum v^ z in the upward (+z)-direction, 
with A a, I a and Ib fixed for a lattice time r of 450 ms. 
The maximum drift velocity clearly increases when Ab 
increases and approaches the F e = 4 state. Continuing 
the increase of Ab close to the F e = 4 state results in an 
abrupt decrease of the TOF-signal due to resonant trap 
losses. However, changing A a, I a and 7b, while keeping 
Ab fixed, affect ip z less significantly and in a more am- 
biguous way compared with changing Z\b, as can be seen 
in figure [8] With the large parameter space and with the 
high degree of coupling between the different parameters, 
it is difficult to extract any unambiguous trends from the 
data. The behaviour of this type of BM is complex and 
since it is working in the dissipative regime, a strong cou- 
pling between friction, heating and dissipation is present. 
Changing any parameter, simultaneously changes the fric- 
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Fig. 7. Maximum vertical drift velocity Vd z as 
of zA B , where A A 33r, I A 4.17 mW/cm 2 , 
mW/cm 2 and the lattice time r = 450 ms. The 
F e = 4 resonance occurs at Ab = 48.T. 
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Fig. 8. Maximum drift velocity Vd z in the (+z)-direction for 
different sets of parameters, plotted as a function of I a- Mean- 
while, Ab = 40.T is kept fixed during a lattice time r = 350 
ms. 



tion, heating and dissipation in the system as well as the 
transfer rates. Basically, since the rectification mechanism 
emanates from the asymmetry in the transfer rates be- 
tween the two OLs, it seems that the asymmetry con- 
ditions for larger induced drifts generally increase when 
approaching the F e = 4 resonance. When detuning closer 
to the F e = 4 transition for lattice B, the possibility for 
optical pumping of atoms to the F g = 3 ground state due 
to absorption of non-resonant photons increases, which 
in turn seem to enlarge the BM effect since atoms are 
spending longer time in the short lived optical lattice A. 
Note that the optimal asymmetry in the transfer rates, 
maximising the induced drift, does not correspond to the 
maximal asymmetry [20] . Typically, the ratio between the 
transfer rates 7a^b : 7b^a are in the order of 9:1, which 
is due to the fact that one is closed and one is an open 
transition [T8] . 




Fig. 9. (Colour online) False colour images of the atomic cloud 
at two different interaction times n and T2, to being with no 
interaction time (r = s) . The final image is acquired after an 
interaction time of t<i ^0.9 s (n ^0.5 s). a) A phase shift only 
along z, showing an upward drift, b) Different phase shift in 
z. A downward drift is evident, ip x and <p y still being zero, c) 
Lp z and Lp y are zero. Only a horizontal drift along x is evident, 
d) Both (p x and ip z are set for a maximum drift in both the z 
and x-directions, resulting in a diagonal drift. 



3.4 Evidence for horizontal and vertical drifts 

Our lattice structures are periodic in three dimensions. 
We can adjust the relative spatial phases (p z , cp y and (p x 
at will. Thus, the Brownian motor works also horizontally, 
and indeed in an arbitrary direction. To confirm this we 
measure the position of the atomic cloud in the xz-plane as 
a function of the lattice time r. This is done by imaging the 
shadow of the atomic cloud, transiently illuminated by a 
weak probe beam, on a CCD-detector. Figure [9] shows false 
colour images of the atomic cloud for different successive 
lattice times. An induced drift is evident both along z 
and x, and as well as in the diagonal xz-dtiection for an 
appropriate choice of <p x and cp z (jp y = 0) which generates 
the largest drift velocities while keeping Z\b, A a, Ib and 
Ia fixed. To confirm the induced drift dependence along 
x as a function of ip x , we determine for each change in 
(f x the centre-of-mass ^-position by a Gaussian fit to the 
images, while cp z and cp y where kept fixed at zero. For 
comparison reason with TOF results in the z-direction, 
the same lattice time r = 350 ms was used for one typical 



set of parameters, see figure 10 



4 Numerical simulations 

In order to understand the qualitative behaviour of our 
BM, we have performed simulations, using a simple model 
displaying the basic characteristics of our system. We con- 
sider a classical Brownian particle which can be in one 
of two internal states, indexed by j, and interacting with 
state-dependent external potentials Uj. In analogy to the 
experimental setup, section [3J we chose Ub to be the light 
shift potential resulting from laser field B (see figure [2]), 
corresponding to the lowest adiabatic optical potential 
[T4] for the F g = 4 — > F e = 5 transition of the 4-beam 
linJ_lin configuration. Consequently, Ua is chosen as the 
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Fig. 10. Experimental results showing the drift velocity Vd x as 
a function of relative spatial phase along x. Here, Ab = 40.T, 
A A 36r, I B 6.13 mW/cm 2 and I A = 6.67 mW/cm 2 
for a lattice time r of 350 ms. The error bars come from the 
uncertainty in the Gaussian fits used to determine the centre- 
of-mass x-position to the images. 




xk,zk 

Fig. 11. (Colour online). One dimensional cuts £7(0,0,2:) 
(blue) and U(x, 0, 0) (red) of the potential surfaces used in 
the simulation, expressed as a function of the light shift |Z\q| 
[13] , where k is the absolute value of the wavevector of the 
laser field. Solid line: lowest adiabatic potential of the F g = 
4 — > F e = 5 transition; dashed line: lowest adiabatic potential 
of the Fv — 3 — > F e — 4 transition. 



lowest adiabatic optical potential for the F g = 3 — > F e = 4 
transition. These potentials are illustrated in figure [TTJ 
where the difference in periodicity and height between the 
x and z directions should be noted (cuts along x or y are 
identical), while Ua and Ub are almost the same. We con- 
sider unequal transfer rates 7a^b 7^ 7b A between the 
two potentials. The Brownian motion is governed by the 
scattering of photons from the optical lattice lasers, rep- 
resented by the diffusion constants D v j(x), while laser 
cooling provides friction. The Fokker-Planck equation for 
the Wigner distribution Wj (x, v, t) of this Brownian parti- 
cle, characterised by its position x, velocity v and internal 




Fig. 12. Results from numerical simulations of the Brownian 
motor mechanism for the vertical z (full line) and horizontal 
x (dashed line) directions. The drift velocity Vd is plotted as 
a function of the relative spatial phase ip = cp A — <£b, varied 
independently along z or x. 



state j as A or B, reads [20] l30] 



[d t +vd x ]W j -d v [v- 



\- VUj(x) + D v j(x)d v ] Wj 
7^(x)W J v- 7j .^,(x)W j (2) 



for j 1 7^ j, with time expressed in units of the inverse of 
the friction coefficient a and position in units of the in- 
verse angular wave vector k of the laser light, such that 
all variables are dimensionless. In addition, we scale the 
potentials by a factor of (45/88) *4.J^] such that the adia- 
batic potential along z, £7(0, 0, z) shown in figure 11 , has a 
depth of w 2.4, allowing a direct comparison with the sine 
potential amplitude A used in ref. [20 . We assume that 
the mean kinetic energy of the system is smaller than the 
depth of the trapping potentials and that the typical fre- 
quency of the oscillations in the potential wells is larger 
than the rate of energy damping, so that the dynamics 
falls into the low-damping regime. Motion is restricted to 
2D (y = 0) and we observe the average velocity in the 
z or x directions as a function of the relative phase shift 
cp in either of these directions. Unless noted otherwise, we 
use arbitrary position-independent values of D v , and 7, 
with 7a^b = 37b^a = 7.5 (i.e., an atom spends 75% of 
its time in the long-lived lattice B). 



Figure 12 shows the dependence of the Brownian mo- 
tor mechanism on the relative phase cp between the two 
potentials, with A = 200 and D v = 75 in both states. 
The overall shape of the curve is similar to what was ob- 
tained for a sine potential in ref. [20] , with a greater drift 
velocity for a phase shift along x, due to the deeper poten- 
tial along that direction. Indeed, figure [13| shows that the 
drift velocity is proportional to the potential depth, up to 
a certain value where the BM effect saturates. Similarly, 
changing the diffusion constant D v affects the BM differ- 
ently depending on the internal state of the atom, as seen 
in figure [14] Increasing the diffusion in both lattices, or 



2 A coefficient of 44/45 appears in the formula for the po- 
tential from the Clebsch-Gordan coefficients for a F g = 4 — > 
F e = 5 transition. 
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Fig. 13. Dependence of the magnitude of drift velocity on 
the depth of the potentials. Circles: identical depth Aa = Ab; 
squares: varying depth of the long-lived lattice Ab for .4a = 
200; triangles: varying depth of the short-lived lattice Aa for 
Ab = 200. (p z = 2tt/3 and D v = 75 in all cases. 




Fig. 14. Dependence of the magnitude of drift velocity on the 
diffusion constant. Circles: identical diffusion D v ,a = ^v,b; 
squares: varying diffusion in the long-lived lattice D v ,b for 
D v ,a = 75; triangles: varying diffusion in the short-lived lattice 
D v ,a for D v ,b = 200. ip z = 2n/3 and A = 200 in all cases. 



only in the long-lived one, results in a decrease of the drift 
velocity, due to the added noise. Conversely, a greater dif- 
fusion in the short-lived lattice slightly increases the drift 
velocity, because the atom can take advantage of the in- 
creased Brownian motion without too much adverse effect 
from isotrope diffusion. 



Compared to the simulations in figure 12, the experi- 
mental results clearly deviate around <p = and 2tt both 
in x and z. This was also noticed in simulations using 
diabatic potential curves [3T] . which prompted us to use 
adiabatic curves in the present study. It is now clear that 
the discrepancy is not due to the shape of the potential, 
which calls for simulations where a more realistic model is 
used to approach the physical system at hand. Moreover, 
we note that the relation between the beam path extension 
and the relative spatial phase, shown in figure [3] indicates 
that Brownian motor effect has the same periodicity as 
the diabatic potential curve [T8] . 



5 Discussion 

There are two ambient forces that may affect the exper- 
imental results. One is due to earth magnetic field, from 
which a slight Zeeman shift could be introduced. This is 
cancelled by a B-field compensation in the experimental 
setup. The second force is due to gravity. It will introduce 
a slight tilt of the lattice potentials, which may results in 
an increased tunnelling probability in one direction. This 
effect is small, and for potential depths used in the ex- 
periment of 100-300£? r . The potential energy L^p = mgh 
contribution due to gravity becomes ~ 10~ 3 E r , where m 
is the atomic mass, g is the gravitational constant and 
h = X/y/2 is the distance between two <r + a - -sites in the 
vertical z-direction (A = 852.3 nm). This clearly indicates 
that this contribution is negligible compared to the poten- 
tial depth. 

While two periodic potentials are used in the model, 
in reality the atoms have magnetic substates, each lead- 
ing to a different light shift potential. A further difficulty 
with our BM is that, since it works in a dissipative regime, 
the parameters, potential depths, magnitude of the diffu- 
sion, and the transition rates between the two lattices all 
depend on the irradiances and the detunings of the OL 
lasers. Changing any parameter results in a different fric- 
tion, diffusion and as well as dissipation. Therefore, one of 
the main features of this paper is to investigate how this 
type of Brownian noise rectifier varies under different con- 
ditions. The complexity of the system makes it difficult to 
directly relate the experimental findings to the simple clas- 
sical model used in the simulation. Quantum mechanical 
simulations in terms of parameters directly controllable 
in experiments, using the full level structure, are under 
development. 

Despite these difficulties, carefully investigations of the 
parameter space, Z\a, Ab-, and /b, allow some general 
conclusions to be drawn. The most striking feature is the 
drastic increase of the effectiveness of the BM when Ab is 
increased. As Ab is increased, the frequency of the light 
will get closer to the F g = 4 — > F e = 4 transition, and 
hence the pumping rate 7b^a increases. Classical simula- 
tions [20 confirm that there is an optimal ratio between 
pumping rates, as long as the inequality remains large. 
Changing Ab will also change the pumping rate between 
the sublevels within the lattice B manifold. Therefore, the 
potential depth, friction and diffusion within lattice B will 
also be modified, both in magnitude and position depen- 
dence, in a non-straightforward way [32] . 

However, in the limit 7b^a/7a^b —> or oo the sys- 
tem is effectively reduced to a single optical lattice, and we 
would expect any BM effect to vanish. Therefore it is clear 
that the optimal BM effect must be achieved for some fi- 
nite ratio 7b^a/7a^b, as has indeed been confirmed by 
classical simulations [20]. This is clearly consistent with 
the data in Fig. [7] Due to experimental difficulties when 
laser B is tuned too close to level A we have not been able 
to extend our data to larger Ab in order to investigate 
if the drift velocity falls off again. Classical simulations 
indicate that this will hapen for 7b^a/7a^b larger than 
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0.1-0.2. Although we do not know the exact relation be- 
tween 7b^a/7a^b and zAb (which will also depend on 
other parameters) measurements of the relative popula- 
tions of the two lattices indicate that this value has not 
been reached [T8] . 

For quantitative analyses, the model using only two 
potentials, while ignoring the internal level structure of 
the two lattices, is not sufficient. This is evident since 
features like the ones at tp = and at cp = 27r for the 
directional variation of the shape of the velocity curve in 
figure [12| do not show up in the experimental graphs (see 
figure [6fT It is also expected that the position dependence 
of pumping, friction and diffusion coefficients will greatly 
affect the shape of the curve. 



6 Conclusion 

In summary, we have demonstrated a Brownian motor 
working in three dimensions with a controllable speed. 
Induced drift velocities in the order of one recoil veloc- 
ity have been achieved. We have showed that some of 
the qualitative features of the experiment can be repro- 
duced theoretically using a classical model similar to that 
in [20] . which has been extended to two dimensions. 

Up to now, a wide variety of BMs have been theo- 
retically investigated and also demonstrated for various 
systems [3j[4] . The main features of our BM can be quali- 
tatively described using a purely classical model shown in 
section [4j Nevertheless, the coupling between the poten- 
tials is driven by quantum jumps as resulting from sponta- 
neous emission, which essentially is a quantum mechanical 
feature. Together with quantised motion, this may open 
the way for creation of a quantum Brownian motor [33] . 
Due to the generality of our scheme, applicability to chem- 
ical and/or biological systems may also be possible. 
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